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Understanding electron behavior in strained graphene 
as a reciprocal space distortion 
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The behavior of electrons in strained graphene is usually described using effective pseudomagnetic 
fields in a Dirac equation. Recently, such description has been completed by the explicit inclusion 
of lattice deformations, although yet it is not clear the physical interpretation of the modified 
resulting equations, or if they really produce new pseudomagnetic fields. Here we consider the 
particular case of a spatially constant strain. Our results indicate that lattice corrections are easily 
understood using a strained reciprocal space, in which the whole energy dispersion is simply shifted 
and deformed by the space shrinking. This leads to a directional dependent Fermi velocity without 
producing pseudomagnetic fields. The corrections due to atomic wavefunction overlap changes tend 
to compensate such effects. Also, the analytical expressions for the shift of the Dirac points as 
well as the corresponding Dirac equation are found. In view of the former results, we discuss the 
range of applicability of the usual approach of considering pseudomagnetic fields in a Dirac equation 
derived from the old Dirac points of the unstrained lattice. Such considerations are important if a 
comparison is desired with experiments or numerical simulations. 

PACS numbers: 73.22.Pr, 81.05.ue 



I. INTRODUCTION 

Since the experimental observation of grapheneji a 
two-dimensional form of carbon, there has been many 
theoretical and experimental studies to understand and 
take advantage of its surprising properties)^"— Among 
its most interesting features, one can cite the peculiar 
interplay between its electronic and mechanic proper- 
ties. Graphene can withstand elastic deformations up 
to 20%, much more than in any other crystal J^ Need- 
less to say, this long interval of elastic response results 
in strong changes of the electronic structure, which of- 
fers a new direction of exploration in electronics: strain 
engineering <i"— The prospect is to explore mechanical 
deformations as a tool for controlling electrical trans- 
port in graphene devices: a technological challenge owing 
to the counterintuitive behavior of electrons as massless 
Dirac fermionsJ^ 

The most popular model proposed in the literature to 
study the concept of strain engineering is based on a com- 
bination of tight-binding (TB) description of the elec- 
trons and linear elasticity theoryi^^rJ^ In this approach, 
where the absence of electron-electron interactions is as- 
sumed, the electronic implications of lattice deformations 
are captured by means of a pseudovector potential A 
which is related to the strain tensor e bj*i^ 



A, 



2a 



{e. 



), Ay^-^{2e^y) (1) 



where a k, 1.42 A is the unstrained carbon-carbon dis- 
tance (see FIG. [T] (a)) and /3 ss 3 modulates the vari- 
ation of hopping energy t of the TB model with the 
changes of the intercarbon distance due to the lattice 
deformations.— 1^ Note that, the x axis is selected parallel 
to the zigzag direction. This /3— dependent pseudovector 



potential gives a coupling of the pseudomagnetic field 
{B = V X A) with the electronic density. The idea of 
pseudomagnetic fields has been key in the understanding 
of the pseudo Landau levels experimental observations 
in strained graphene, which had been theoretically pre- 
dicted earlier J^iii 

In recent works, the standard description of the strain- 
induced vector field has been supplemented with the ex- 
plicit inclusion of the local deformation of the lattice 
vectorsjiSii^ After accounting for the actual atomic po- 
sitions to the TB Hamiltonian, Kitt et al. proposed 
an extra pseudovector potential which is /3— independent 
and different at each of the strained Dirac points. ^^ The 
possible physical relevance of the extra /?— independent 
term, predicted in the work of Kitt et al., as discussed 
by de Juan et al. within the TB approachj^^ They 
also obtained a extra /3— independent pseudovector po- 
tential but with zero curl. Therefore, they concluded that 
in strained graphene no /3— independent pseudomagnetic 
field exist ,^ 

In this paper, we analyze the most simple case, a spa- 
tially uniform strain. The reason is that such case must 
be contained as a limiting case in any of the general theo- 
ries, and at the same time, as shown here, it can be solved 
exactly. Thus, it is an important benchmark tool to com- 
pare with previous approaches. For example, this leads to 
a simple explanation for the lattice correction terms and 
direction dependent Fermi velocity, since both are due to 
the effects of strain in reciprocal space. Hopefully, this 
will help to derive the consequences of lattice corrections 
of flexural modes or curved graphene. The layout of this 
work is the following. In Section II we present the model 
and find the corresponding energy dispersion surface. In 
section III, we discuss the properties of the energy dis- 
persion and find the analytical expressions for the shift 
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FIG. 1. (Color online) (a) Unstrained graphene lattice show- 
ing the vectors Si that point to the neighbors of type A sites, 
(b) the same lattice under a uniform stress, and (c), the 
first Brillouin zone of the reciprocal lattice for the unstrained 
(dashed lines) and strained (solid lines) graphene. Notice how 
the reciprocal lattice is contracted in the direction where the 
lattice is stretched, and the change of the Ko symmetry point 
into K. In (d), it is shown how the distortion of the reciprocal 
lattice transform the original Dirac cone (left) into a distorted 
one (right) with a directional dependent Fermi velocity. 



of the Dirac points and the strained Dirac Hamiltonian. 
Section IV deals with the problem of how the usual pseu- 
domagnetic fields approach needs to be added with some 
requirements in order to compare with experiments and 
simulations. Finally, in the last section, the conclusions 
are given. 



II. MODEL: ELASTICITY & TIGHT-BINDING 

We are interested in uniform planar strain situations, 
i.e., the components of two-dimensional strain tensor e 
are assumed to be position- independent. In this case, the 
displacement vector u{x) is given by u{x) = e ■ x, and 
therefore, the actual position of an atom x' — x + u{x) 
can be written as x' = {I-\-e)-x, being / the 2x2 identity 
matrix. In general, if r represents a general vector in the 
unstrained graphene lattice, its strained counterpart is 
given by the relationship r' = (/ -f e) • r. 

We investigate the electronic implications of strain by 
means of the nearest-neighbor tight-binding Hamiltonian 



if = - ^ ix',na^/fex'+5; + H.c. 



(2) 



where x' runs over all sites of the deformed A sublat- 
tice and 8'^ are the three nearest neighbor vectors. The 



operators a^, and hx'+S' correspond to creating and an- 
nihilating electrons on the sublattices A and B, at sites 
x' and x' + 8'^, respectively. The dispersion relation arise 
upon writing Eq. ^ in the momentum space. For this 
purpose, we replace the creation/annihilation operators 
by their Fourier expansions^! 
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h^.+S'^ = ^ E e-''^-(-+^"+"(-+''"))6fc„ (3b) 

^ fc2 

where N is the number of elementary cells. In Eq. ([2]) 
we have written the hopping integral tx\n as position- 
dependent, but in the considered case (uniform strain), 
it does not depend on the position, only on the direction: 

Under these considerations, the calculation of the 
Hamiltonian in the fc-space is fairly straightforward, H 
becomes 



H 






t„e-''=(^+^)-^"ajt6fc + H.c. 



(4) 



From this equation, it follows that the dispersion relation 
of graphene under spatially uniform strain is 



E(k) 



±|^t„e-^'=-(^+^)"'" 



(5) 



Notice that in the general case of a non-uniform strain, 
u{x) is not given as u{x) = e{x)-x instead e{x) needs to 
be replaced by the displacement gradient tensor Vii J^ 



III. ENERGY SPECTRUM OF STRAINED 
GRAPHENE 

The variation of hopping energy t„ with the changes 
of the intcrcarbon distance fulfills a physically accurate 
exponential decay i„ = to exp[— /3(|5^|/a — 1)], with 
to ~ 2.7 eV being the equilibrium hopping energyj^i^ 
Nevertheless, for the sake of comparison with previous 
works we consider first order in strain. 



tr, 
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Defining the three nearest neighbor vectors as depicted 
in FIG. n 

^i = ^(V3,l), 52-^(V3,l), <53 = a(0,-l), (7) 

and plugging Eq. (jH]) into Eq. ([5]), we get the following 
expression for the dispersion. 



Eik) = ±to V3 + /(fc) - /3 [3tre + /,(fc)] + /32/e^(fe), 



(8) 



where f{k) has exactly the same functional form of its 
counterpart in unstrained graphenepi 

/(fe)^2cos(\/3fc»+4cos(^-^^)cos(^), (9) 

but now evaluated in different points of reciprocal space, 
since here k* — (fc*, k*) is given by the transformation, 



k* = {I + e)-k. 



(10) 



This last equation is very important. It provides a 
mapping of the original reciprocal space into a new dis- 
torted one. As we will see, this mapping and the fact 
that /(fe) is equal to its undistorted counterpart, leads 
to pure geometrical effects that only very recently have 
been identified. The other terms depend on the same 
distortion, but contain hopping corrections. These terms 
are explicitly detailed in the appendix, fdk) contains 
the modification of the spectrum due to first order in /3, 
while /g2 (fc) is the second order correction in /3. 

Several important consequences follow from these 
equations. First of all, we observe that in the case of 
deforming the lattice without changing the hopping pa- 
rameters, i.e., if we deform the lattice keeping /3 = 0, we 
get that. 



Eik)^±h^S + f{k). 



(11) 



This corresponds to the same dispersion relationship 
observed in graphene, but now with different reciprocal 
vectors, which are obtained by applying strain to the 
original reciprocal vectors. In other words, the space is 
strained while the eigenvalues remain the same. As a 
consequence, the Dirac cone just changes its shape due 
to this lattice deformation, as illustrated in FIG.[T] This 
is exactly the result that we would obtain if a diagonaliza- 
tion of the tight-binding Hamiltonian is performed using 
a computer. Since /3 = and the connectivity of the lat- 
tice is not changed, the eigenvalues of the Hamiltonian 
must remain equal to the undistorted case. Only when 
a plot is made against the wavevectors, the cone turns 
out to be distorted, as shown in FIG.[T](c) and (d). For 
example, the brick wall lattice has the same connectivity 
as graphene, and thus the spectrum must be the same. 
However, only when the spectrum is plotted in recip- 
rocal space, the energy-momentum relationship appears 
distorted. 

Also, the case (3 — allows to appreciate a subtle 
point. Since the spectrum is the same than in unstrained 
graphene, it is easy to see that the K symmetry points 
of the distorted lattice coincide with the Dirac point of 
the new relationship given by Eq. (jll[) . In other words, 
the condition E{Kd) — 0, which defines the Kd Dirac 
points, corresponds to Ko — K , where K is the image 
of the point Kq under the mapping K = {I + e)^^ ■ Kq. 
Thus, for /3 = it makes sense to develop the TB Hamil- 
tonian around the original Dirac points, as fe = K^ + q, 
with |q| ^ \Kd\- It is easy to show that the pure geo- 
metrical distortion allows to write the Dirac Hamiltonian 



H = vqct ■ {I + e) ■ q = vqct' ■ q 



(12) 



being q the momentum measured relatively to the Dirac 
points, Wo = 3tQa/2 the Fermi velocity for the unde- 
formed lattice, cr = (axTTy) the two Pauli matrices and 
cr' = {I + e)-a turns out to be the spinorial connectiouj^ 
From this equation follows a directional dependent Fermi 
velocity for strained graphene, which has also been found 
in other worksf^^ Furthermore, a directional dependent 
velocity appears simply by looking at the isoenergetic 
curves of Eq. pT|) around K. In this case one obtains 



Eiq)^ =. [voa- •(/ + £)• qf 



(13) 



where q is the nromentum measured relatively to the 
Dirac point, therefore, the isoenergetic curves around K 
are rotated ellipses, as depicted FIG.[T](c). 

At this point, we can conclude that the basic mech- 
anism behind the anisotropic Fermi velocity is the dis- 
tortion of the reciprocal space. Such distortion gives a 
simple interpretation to the resulting geometric crystal 
frame terms that appears in the covariant version of the 
equations.^- Clearly, there are not associated pseudomag- 
netic fields^iiS. 

Let us now consider the case in which the space is dis- 
torted and the hopping is changed, i.e., /? 7^ 0. Here, we 
will have two effects. Again we have the pure geometri- 
cal distortion due to the strain of the reciprocal space, 
but at the same time, there is a change in the spectrum. 
This last effect is the only one observed when a diago- 
nalization of the Hamiltonian is performed in a computer 
for a finite number of atoms. In FIG. [5] , a comparison 
between the cases /3 = and /3 7^ is presented for E{k). 
FIG. [21 was made for a zigzag uniaxial strain of 5%, i.e., 
^xx — 0.05, e^y — and tyy = —crexx, being a the Poisson 
ratio, which is very small for graphene, a ^ 0.1 — 0.15, 
according to some theoretical estimations. ^^'^'^ 

As can be seen, the effect of /? 7^ is to distort the 
/3 = case in such a way that it tends to compensate the 
strain of the reciprocal space, i.e., the ellipses are rotated 
by 7r/2 for a realistic value of /3. The physical reason for 
this to happen is that a stretched direction in real space 
shirks in reciprocal space, resulting in a higher Fermi 
velocity, while in the same direction, the orbital overlap 
decreases since the distance between atoms increases (see 
FIG. [1]). This tends to reduce the Fermi velocity. As 
a result, lattice distortion and hopping changes tend to 
compensate. This fact can also be seen in the movement 
of the Dirac points. From the FIG.Hl one can see that the 
Dirac points for the case /3 7^ are closer to the original 
ones than their (3 — counterparts. 

The qualitative results discussed above, and depicted 
in FIG. m can be understood by finding analytical ex- 
pressions for K]j and H. The position of Ko can be 
obtained from the condition E{Ko) = 0. Up to first 
order in strain, Kd is given by, 

KD^{I + e)-^-{Ko+^A)c^K + iA (14) 





2.3032 

FIG. 2. (Color online) Isoenergetic curves obtained from the energy dispersion in reciprocal space obtained from Eq. ((8]). A 
blow up is presented around the Dirac points Kd for each surface. Case (a) corresponds to unstrained graphene, (b) to strained 
graphene with e^cx ~ 0.05, e^y = 0,eyy — —atxx, and (c) to strained graphene with P k, Z, and the same strain tensor than 



in case (b). Notice how although the strain tensor is the same in cases (b) and (c) 
reciprocal space deformation and hooping effects tend to compensate. 



the ellipses are rotated by 7r/2, since the 



with A given by Eq. ([T]) and ^ the valley index of 
Kq^ The previous equation confirms the remark that 
the Dirac points for /3 ^ do not coincide with the K 
high symmetry points of the strained Brillouin zone. The 
shift, which is only produced by /3, is given by the pseu- 
dovector potential, and do not depend on Kq- 

Furthermore, once the points Kjj are known, it is pos- 
sible to obtain a new Dirac Hamiltonian considering the 
lattice correction and orbital overlap changes. To do this, 
one can develop Eq. ([5]) around the new Dirac points us- 
ing Eq. (|14p. and derive that 

H = vo(T-{I + e-Pe)-q (15) 

which is a general version of Eq. p^ . since /3 effects are 
included. Notice that the isoenergetic curves around Kd 
remain being ellipses, as depicted in FIG.[2](c), but with 
different values of the semi-axes owing to the /3 correc- 
tions. The last equation clearly shows the tendency of (3 
to cancel the lattice corrections. Eqs. (1141) and (flS)) are 
among the main contributions of this paper. 



IV. EXPERIMENTAL OBSERVATION OF 
PSEUDOMAGNETIC FIELDS 

Prom the point of view developed in the previous sec- 
tion, is clear that basically, the Dirac cone is translated 



and distorted. As a result, if one tries to derive an effec- 
tive Dirac equation using Kq as starting points to develop 
E{k) as k — Ko + q, the resulting energy can be quite far 
away from the Fermi energy, as shown in FIG. [3] This 
poses a problem that has been overlooked in the usual 
treatment of strain in graphene using pseudomagnetic 
fields in the Dirac equations. 

In general, if Eq. ^ is developed around a general 
point in reciprocal space given by Kg, we get, 

E^{KG+q) « E\KG)+VE\KGyq+y-VVE^iKGyq 

(16) 
where 'VE'^{Kg) is the Jacobian vector and 
WE'^{Kg), the Hessian matrix of E'^{k), which 
are evaluated at A; = Kg- 

In the usual procedure Kg — Kq. However, E^{Kf)) ^ 
and 'VE^{Kg) y^ 0. This produces an energy shift 
and a q dependent term, observed in other approaches^, 
which complicates somehow the description of the dy- 
namics. 

This also poses an issue concerning the experimental 
possibility of observing the pseudomagnetic fields. Since 
the energy evaluated at the original Dirac point E{Ko) is 
different form zero, the Fermi energy does not fall in this 
point, as we illustrate in FIG. |31 In general, if a exper- 
iment is performed at temperature T, and the chemical 




FIG. 3. (Color online) The Dirac cone in unstrained (dashed 
line) and strained graphene (solid line) and the experimental 
observation of electron behavior for a probe that shifts the 
chemical potential (n) with respect to the Fermi energy. The 
shaded box indicates the width of the thermal selector due 
to the Fermi-Dirac distribution. The effective Dirac equation 
with pseudomagnetic fields can be obtained by developing 
around the original Kq points, or in the Dirac points Kd of 
the strained lattice. For n — 0, only the last approach will 
work for low temperatures. 



potential /i is shifted by a field, the condition to observe 
the pseudomagnetic fields in the usual derivation around 
the original Dirac point must satisfy. 



E{Ko) - fi \< keT 



(17) 



since the difference between E{Kq) and /i must be less 
than a zone defined from the derivative of the Fermi- 
Dirac distribution against the energy, as explained in 
FIG. [3] using a box around ^. As T — ?► 0, the deriva- 
tive is a delta function centered around the Fermi en- 
ergy, and the pseudomagnetic fields calculated from Kq 
are usually far from the region of validity. For exam- 
ple, even a zigzag uniaxial strain of 1%, will produce a 
E{Kq) > 27 meV, which is much higher than the ther- 
mal width of ksT sa 8.6 meV, obtained at T = 10 K. 
This breaks the approximation of using pseudomagnetic 
fields in a Dirac equation unless a very well defined field 
is used. 

The option is to have a better description of the en- 
ergy dispersion near the Fermi energy, by developing Eq. 
([T6| around the true Dirac points of the strained lattice, 
i.e., by setting Kg = K^, for which the corresponding 
energies fall at the Fermi level. In this case, 

E''{KD + q)~\q-VVE''{KD)-q (18) 



since E'^{Kd) = and VE'^{Kd) = 0. Now one obtains 
an energy dispersion which corresponds to a distorted 
cone, with a directional dependent Fermi velocity given 
by the elements of the Hessian oiE^{k) evaluated at Kn . 
This result is the same than the one obtained from the 
Dirac Hamiltonian given by Eq. (|15|) . 



V. CONCLUSIONS 



In conclusion, we have analyzed the case of a spatially 
uniform strain in graphene. The lattice correction terms 
are simply an effect of a strained reciprocal space. As 
a consequence, the Dirac cones are deformed and trans- 
lated. No pseudomagnetic fields are associated to such 
terms, as has been recently discussed^ii^. When hopping 
changes are considered, there is an extra deformation of 
the cone that tends to cancel the effect of the reciprocal 
space strain. The new Dirac points of the strained Hamil- 
tonian do not coincide with the K symmetry points of 
the strained reciprocal lattice. Due to this fact, the effec- 
tive Dirac equation can be obtained either by developing 
around the old or new Dirac points. If the old points 
are chosen, as usual in graphene's literature, there is a 
restriction to observe the dynamics produced by the cal- 
culated pseudomagnetic fields since only for very high 
temperatures or carefully designed probes it is possible 
to compare with the usual theory. In computer simula- 
tions, it is also important to distinguish between lattice 
distortion effects and connectivity matrix. Some of these 
issues, can explain differences between theory and simu- 
lations in graphenei^ 

Finally, it is worthwhile mentioning that although we 
only treated a particular case, the ideas and lessons ob- 
tained from this study can be translated to general cases, 
as we will show in forthcoming works. 
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Appendix 

In this appendix, we provide explicit expressions for 
the last terms in Eq. ([8]), 



fe{k) = {Sexx + e.yy) cos(%/3fc*a) + [ie^x + ^^yy) cos( — -^) cos(— ;|-) - 2V3ea;y sin( 



vSfcJa^ 



3fc,*a 



,V3fc* 



I sin(- 



Zkla 



U{k) = -(9eL + Sexxey,, + 96^ + Uel + ((3e,, + eyyf 



yy v"^^^ 



I cos 



12e^j^ ) cos 



(VSfc^) 



(%/3fc*/2)cos(3fc;/2) - 16V3eyj,e:,j;Sin(V3fc;/2)sin(3fc;/2)). 
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